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Abstract 

We present here some striking discrepancies in the results of ac and dc magnetization measurements performed in a single 
crystal of low Tc superconductor, Yb 3 Rh 4 Sni 3 . Fingerprint of a transition from an ordered vortex lattice a la Bragg glass (BG) 
phase to a partially-disordered vortex glass (VG) like phase gets unearthed under the influence of an ac driving force present 
inevitably in the isothermal ac susceptibility measurements. In contrast to its well-known effect of improving the 

state of spatial order in the vortex matter, the ac drive is surprisingly found to promote disorder by assisting the BG to VG 
transition to occur at a lower field value in this compound. On the other hand, the isothermal dc magnetization (M-H) scans, 
devoid of such a driving force, do not reveal this transition; they instead yield signature of another order-disorder transition at 
elevated fields, viz., peak effect (PE), located substantially above the BG to VG transition observed in x'(H) runs. Further, 
the evolution of PE feature with increasing field as observed in isofield ac susceptibility {x'{T)) plots indicates emergence of an 
ordered vortex configuration (BG) from a disordered phase above a certain field, H* (~ 4kOe). Below H*, the vortex matter 
created via field - cooling (FG) is found to be better spatially ordered than that prepared in zero field - cooled (ZFC) mode. This 
is contrary to the usual behavior anticipated near the high-held order-disorder transition (PE) wherein a EC state is supposed 
to be a supercooled disordered phase and the ZFG state is comparatively better ordered. 

PACS numbers: 74.25.0p, 74.25.Ha, 74.25.Dw 
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I. INTRODUCTION 

In the context of the mixed state of a type - II super¬ 
conductor, the seminal discovery [Il,i of a well-ordered 
thermodynamic phase, viz., Bragg glass (BG) exhibiting 
Bragg’s reflections, and its possible transition(s) to a 
disordered phase devoid of Bragg’s reflections had lead to 
a generic vortex phase diagram applicable to almost 
all pinned superconductors. Owing to the possibility of 
a sudden proliferation of dislocations on progressively in¬ 
creasing the magnetic field at a constant temperature, the 
quasi - long range ordered BG phase is anticipated [l|, 0 
to transform first into a multi - domain (partially disor¬ 
dered) vortex glass (VG) phase. Such a transition usually 
reflects as a second peak in (isothermal) magnetization 
{M-H) loops, termed as the second magnetization peak 
(SMP) anomaly [iOl. Thereafter, at elevated fields 
closer to the upper critical field {IIc 2 {T)), there occurs 
another anomaly, known as the quintessential peak effect 
(PE) phenomenon SiEIIml in field/temperature 
variation in critical current density, jc{II,T). The PE is 
argued [l^ to signal the collapse of the elasticity of an 
ordered vortex lattice at a rate faster than the pinning 
force density near IIc 2 {T). 

Although theoretical treatment related to the BG to 
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VG transition exists in the literature (for example, as in 
Ref. i). however, the experimental tools, particularly 
those employed to explore the bulk pinning properties, 
such as, the dc magnetization (M-H) measurements do 
not always capture this transition. Therefore, in the con¬ 
text of vortex phase diagram studies, it is tempting to ask 
a question- Is the BC to VC phase transition generic? If 
it is so, then this anomaly must get exposed in the H-T 
space of all pinned superconductors, possessing a cer¬ 
tain amount of quenched disorder. To address this issue, 
we have investigated via detailed magnetization measure¬ 
ments, a low Tc superconductor Yb 3 Rh 4 Sni 3 , which has 
so far been reported to display only the PE phe¬ 

nomenon. A claim of the presence of SMP anomaly (or 
its counterpart) in this compound has so far not been 
made by anyone. We have now found that the BG to VG 
(akin to SMP anomaly) transition in this compound gets 
exposed prior to the onset of PE under the influence of an 
ac driving force present in the (isothermal) ac suscepti¬ 
bility scans {x'{H)). Gounter-intuitively, shaking of the 
vortex array by an ac drive in the present study has been 
seen to promote the spatial disordering in the vortex mat¬ 
ter by assisting the BG to VG transition process. This 
observation is in complete contrast to the usual role of an 
ac driving force, which is to improve the state of spatial 
order in multi - domain vortex matter as reflected [l7j| by 
the enhanced brightness of the Bragg spots in the field- 
temperature phase space prior to crossover to the PE 
region in the vortex phase diagram. Another interesting 
aspect of the present study is a revelation of the inequal¬ 
ity {H) < at lower fields (below a char¬ 

acteristic field, H* 4kOe) . This amounts to stating 
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that the vortex matter created in the field - cooled mode 
(FC) below H* exhibits better spatial ordering than that 
obtained in the zero field-cooled manner. This feature 
also is in sharp contrast to the behavior seen at 

the higher fields, where one encounters order - disorder 
transition (a la PE phenomenon). The vortex matter 
created via field - cooling is expected to be spatially more 
disordered (than the ZFC state) due to supercooling of 
an amorphous vortex matter below the PE region. We 
present evidences that H* signifies a crossover regime 
from an ordered BG phase into a disordered amorphous 
phase while reducing the field. 

II. EXPERIMENTAL DETAILS 

Single crystals of Yb 3 Rh 4 Sni 3 have been grown by tin 
flux method [1^. The specimen chosen for the present 
study is platelet shaped, with a planar area of 5.76 mm^ 
and thickness of 0.62 mm. The superconducting tran¬ 
sition temperature (Tc) of this crystal is found to be 
nearly 7.5 K. Magnetization data, both ac as well as 
dc, were recorded using the same instrument, viz., Su¬ 
perconducting Quantum Interference Device-Vibrating 
Sample Magnetometer (SQUID-VSM, Quantum Design 
Inc., USA). The magnetic field was directed along the 
crystalline [110] axis, with a possible error in the align¬ 
ment to be within 5^^. The demagnetization factor in 
this orientation is expected to be small as the field is 
applied in the plane of the thin platelet, i.e., normal to 
the smallest dimension (thickness) of the sample. In the 
dc measurements, we kept the amplitude of vibration of 
the sample to be small (~ 0.5 mm) so as to minimize the 
field-inhomogeneity along the scan length. During the 
ac susceptibility measurements T)), an ac field of 

amplitude 1 Oe and frequency 211 Hz was superimposed 
on the applied dc field. 

III. RESULTS 

A. Isothermal dc M—H loops: Manifestation of 
the peak effect phenomenon 

The inset panel of Fig. 1(a) displays the first two quad¬ 
rants of an isothermal dc magnetization hysteresis loop 
(M-H) obtained at T = 2K for field applied parallel 
to the [110] plane of Yb 3 Rh 4 Sni 3 crystal. The M{H) 
curve can be seen to be hysteretic between the forward 
{M{H^)) and the reverse {M{H~)) sweeps of the mag¬ 
netic field, as expected for a pinned type - II superconduc¬ 
tor. However, there exists an unusual enhancement in the 
hysteresis width (AM(i7) = M(i7+) — M{H~)) prior to 
Hc 2 , as apparent from the encircled portion of the M{H) 
curve. A magnified view of this portion is displayed on 
an expanded scale in the main panel of Fig. 1(a). Using 
a prescription of the Bean’s Critical State Model [2^, Fi- 
etz and Webb have shown [2^ that the hysteresis width 



Figure 1: (Color online) (a) An expanded portion of (isother¬ 
mal) dc M-H loop at T = 2 K demonstrating a typical PE 
phenomenon in a single crystal of Yb 3 Rh 4 Sni 3 . The inset 
displays the first two quadrants of this loop in the entire field 
range of investigation, i.e., 0 < H < 28kOe. (b) The M-H 
curves at different temperatures show the evolution of the PE 
feature as it occurs at different fields. 


(AM) can be taken as a measure of the critical current 
density, jc{H,T). Therefore, the enhancement in AM 
reflects an unusual increase in jc{H) a little below Hc 2 , 
which can be identified as the PE phenomenon. The on¬ 
set field and the peak field (H^e^’'*") of the PE 

stand located in the main panel of Fig. 1(a). The merger 
of M(i7+) and M{H~) beyond the bubble feature iden¬ 
tifies the irreversibility field (Hirr)- 

We show in Fig. 1(b), the fingerprint of PE feature 
at different temperatures, as indicated. On increasing 
the temperature, a systematic decrease in the onset and 
the peak field values of the PE is well apparent. The 
decrease in and with T, nearly following 

the variation of Hc 2 with T, is a characteristic feature 
of the PE phenomenon. The order-disorder transition 
pertaining to PE has been argued to have first - order 
character, as was well evident in the results of small angle 
neutron scattering study in a low Tc superconductor Nb 
[13 and those of scanning Hall probe microscopy in 2H- 
NbSea 
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B. Isothermal ac susceptibility x'{H) responses: 
Identification of an additional anomaly in conjunction 
with the PE phenomenon 

To investigate further the order - disorder transition(s) 
in the vortex matter, we recorded the isothermal ac sus¬ 
ceptibility {x'{II)) responses as shown in Fig. 2. The 
sample was initially cooled - down to a chosen tempera¬ 
ture in (near) zero field and thereafter the x' data were 
recorded while ramping the field to higher values. The 
lx'I values at T = 2K can be seen to fall monotonically 
with increasing field until about a certain field, marked as 
^on,ac 12 5kOe, cf. main panel of Fig. 2(a)). Above 
Hon,ac^ a broad (anomalous) dip-like feature encompass¬ 
ing a large field interval (i.e., < H < Hc 2 ) can 

be noticed. Note that the identification of onset field 
(^Hon,ac) iijjg jjp feature was made possible by locating 
the first zero-crossing in the derivative plot of x'{H) as 



H (kOe) 

Figure 2: (Color online) (a) The real part of (isothermal) ac 
susceptibility (x'iH)) data plotted against the dc magnetic 
field at T = 2K. A broad dip-like anomaly triggering at a 
field marked as has been observed. The inset panel 

illustrates a portion of dy'/dH curve at T = 2K which helps 
to identify the onset (JFp"’“‘^) and peak (/Ip^’“°) positions of 
the anomaly identified respectively via the first and second 
zero - crossings of the derivative plot, (b) x'{H) curves at 
different temperatures show how the anomaly changes from a 
broad to a narrow one with increase in temperature. 


displayed in an inset panel of Fig. 2(a). The second zero- 
crossing in dy'/dH identifies the peak position (marked 
as of the anomalous variation in x'{II) across 

the dip region. We recall here that the x ' value reflects 
jc through the two relations [l^l, (i) x^ + cthac/jc 

for hac < h* and (ii) y' ~ -pjc/hac if hac > h*, where 
a and /3 are size and geometry dependent factor, hac is 
the amplitude of the ac field and h* is the ac field for 
full penetration. Following equation (i), the anomalous 
variation in x'{II) above in Fig. 2(a) amounts to 

an unusual increase in otherwise monotonically decreas¬ 
ing jc{H). Similar anomalous behaviour in x'{H) can 
also be observed at different temperatures as illustrated 
in Fig. 2(b). While the dip feature remains quite broad at 
lower temperatures, one can note that the width of this 
anomalous region gets substantially reduced on higher 
temperature side. For example, the dip in x'{H) at 5K 
is observed to be much sharper with a narrow transition 
width (of nearly 2 kOe) prior to approaching Hc 2 which, 
following equation (ii), amounts to a rapid increase in 
jc{H) echoing the characteristics of the quintessential PE 
phenomenon. 

A comparison of x'{H) data presented in Fig. 2(a) 
with that of the M-H loop shown in Fig. 1(a), tells 
us that the peak field value rs 20kOe) of the 

anomaly observed in the former is nearly the same as the 
corresponding peak field value (Ff^g’^*^) of the PE ob¬ 
served at the same temperature in the latter. However, 
if we seek the analog of onset field of the anomaly seen 
in ac x'(-ff), viz., (« 12.5kOe, Fig. 2(a)), we find 

that no anomalous feature stands depicted at this field 
value in the dc M-H loop (cf. Fig. 1(a)). This implies an 
additional anomaly has got unearthed (at iF°"’“°) in the 
x'{H) scans in the specimen of Yb 3 Rh 4 Sni 3 which is not 
observed in its M-H runs. To seek the rationalization 
of this discrepancy in the ac and dc magnetization data, 
one can consider the presence of an ac driving force in¬ 
evitably superimposed on the dc field during the x'{H) 
runs as an important difference between the two sets of 
measurements. Two important factors in conjunction, 
viz., (i) the shaking effect of an ac driving force on the 
vortex matter and, (ii) driving effect due to a continu¬ 
ous ramping of the dc magnetic field, may be acting as a 
combined driving force leading to the triggering of addi¬ 
tional anomaly in x'{H) runs. To shed some more light 
onto this assertion, it became tempting to explore the 
isofield ac susceptibility responses, where an effect due 
to the ramping of magnetic field would not be present. 
These are described ahead. 


C. Isofleld ac susceptibility X (T) responses: iden¬ 
tification of the PE phenomenon 

Figure 3 displays the temperature dependences of the 
in-phase ac susceptibility (y'iT)) obtained at various 
fixed dc fields. The sample was initially cooled down to 
1.8 K in (nominal) zero field, a desired field was then ap- 
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Figure 3: (Color online) Observation of PE anomaly in tem¬ 
perature-dependent ac susceptibility {x'{T)) runs at various 
fixed dc fields. Onset and peak (Tp‘^) temperature of 

PE have been marked in the x'{T) curve for Hdc = ISkOe. 
The inset shows the occurrence of zero-field superconduct¬ 
ing-normal transition at T « 7.5 K. 


plied and the x'(T) data were recorded while warming up 
to higher temperatures (T > Tc). The zero-field super¬ 
conducting transition temperature, TcCO) identified via 
the onset of diamagnetic response is found to be about 
7.5K (cf. inset panel of Fig. 3). Fingerprint of the PE 
feature identified by a dip - like characteristic can be ob¬ 
served in all the curves shown in the main panel of Fig. 3. 
At H = 5kOe, the dip in x'(r) is found to be less 
prominent, however, it evolves into a sharp (negative) 
peak at higher field values (see, e.g., x'(T) response in 
H = lOkOe). On further increasing the field to ISkOe, 
the PE region becomes broad, which is apparent from 
a wider gap between the onset (r°”’“°) and the peak 
(Tp'^) temperatures of the PE marked in the main panel 
of Fig. 3. Such a broadening seen in the PE feature at 
higher fields may be ascribed to effects of an enhance¬ 
ment in ’effective pinning’ with field, as articulated by 
Giamarchi and Le Doussal 0,1 and elucidated in a crys¬ 
tal of low Tc superconductor 2H-NbSe2 by Banerjee et 

ai, [m. 

A lesser developed (PE) anomaly at lower field value 
(~ 5 kOe) suggests the possibility of a lesser ordered vor¬ 
tex matter emerging at that end as well. Such a trend is 
consistent with the notion that an ordered vortex lattice 
(Bragg glass) undergoes a transition into a disordered 
phase on lowering the field [H, [2^. 

Across the field interval from about 5kOe to 14kOe, 
the vortex matter seems to evolve into a better spatially 
ordered (Bragg glass) phase, as is evidenced by a sharp 
PE feature at lOkOe in Fig. 3. 



Figure 4: (Color online) A sketch of vortex phase diagram in 
our crystal of Yb 3 Rh 4 Sni 3 . H°”’“‘=(r) and have 

been extracted from x'{H) and M-H plots, respectively. Be¬ 
low T = 4.5 K, the former resembles the BG to VG transition 
line (refer text in section-4, “Discussion”) while the latter 
portrays the characteristics of onset of the PE anomaly. In 
the boxed - region, both J7°"’““(T) and //“"’“^“(T) lines be¬ 
have as the onset of the PE transition. The onset temper¬ 
atures of the PE obtained from isofield x'(T) 

scans almost fall on the onset field values of the 

PE extracted from M-H loops. Hc 2 {T) taken from M-H 
loops depicts usual linear fall with increase in temperature. 


D. H—T phase diagram 

It is instructive to compare the results of ac and 
dc magnetization data of Figs. 1 to 3. For this pur¬ 
pose, we present in Fig. 4, a H-T phase diagram of 
Yb 3 Rh 4 Sni 3 , which comprises the field/temperature val¬ 
ues corresponding to the onset positions of the anomalies 
seen in Figs. 1 to 3. The following features in this phase 
diagram are noteworthy: 

1. The onset field values (i/°”’‘^^(T)) of the PE (open 
triangles) acquired from the M-H data (Fig. 1) fall 
smoothly with the increase in temperature, a trend sim¬ 
ilar to the Hc 2 {T) line. 

2. The corresponding onset field values (Ff°”’“'^(T) 

shown by closed triangles) of the anomaly observed in 
x'{H) plots (Fig. 2), however, do not coincide with 
the onset position (iJ™’‘^'^(T)) of the PE over a signif¬ 
icantly large portion of the H-T space (i.e., T < 4.5 K). 
Here, the i7°"'’“°(T) line stays substantially below the 
Hon,dc^rj.) Moreover, iJ°"-“^(r) values exhibit a 

weaker temperature dependence at temperatures below 
4.5 K unlike the faster temperature dependence seen in 
case of H°ff‘^ffT). The two onset field lines, H°ff‘^%T) 












D H-T phase diagram 


5 


and come closer to each other only in the 

boxed region (i.e., at T > 4.5 K and 5 kOe < H < 10 kOe) 
and show identical temperature dependence. 

3. The discrepancy in the phase diagram (Fig. 4) per¬ 
taining to the location of onset helds of the anomalies 
seen in the (isothermal) M-H and x'{H) data prompts 
the need to take into consideration the results of isoheld 
x'(r) scans (Fig. 3) as well. For this purpose, we have 
plotted in the phase diagram the onset temperatures 
( 2 ^n,oc(^)) pg (shown by stars) obtained from 

x'{T) scans. It is curious to note that the T°"’“'^(iJ) val¬ 
ues fall almost on the onset field (iJ°”’‘^°(T)) line of the 
PE transition. We draw an important inference here that 
the results of isothermal dc M-H loops and the isofield 
ac x'{T) scans predict the onset of the PE anomaly at 
nearly the same phase boundary, i.e., 
line. 

Clearly, there exists an additional anomaly (apart from 
the PE) located deeper (at iJ°"’’“'^(T)) in the mixed state 
of Yb 3 Rh 4 Sni 3 , which has got unveiled only from the 
outcomes of the x'{H) data. The other two measure¬ 
ment techniques (viz., the M-H and x'{T) scans), on 
the other hand, reveal the fingerprints of only one kind 
of anomaly, i.e., the PE. This corroborates our previous 
assertion that both an ac driving force as well as the 
dc magnetic field ramping together trigger the said ad¬ 
ditional anomaly well below the onset of PE (these two 
factors not being present together during the individual 
dc M-H and ac x'(T) runs). 

The lowest field down to which the PE could be dis¬ 
cernible in Fig. 4 is nearly 4.5 kOe. Since there is hardly 
any fingerprint of PE feature below this value, we could 
surmise that the vortex matter there may be disordered. 
This proposition is explored further via the thermomag- 
netic history - dependent magnetization measurements. 


E. Thermomagnetic history dependence in 
jc{H,T): FC state is more ordered than the ZFC at 
low fields 

Figure 5 displays the shielding responses {x'{H)) ob¬ 
tained at 2 K for the system prepared in two different 
histories, viz., zero field-cooled (ZFC) and field-cooled 
(FC) modes. The isothermal x'{H) data in the ZFC 
mode (open circles) are the same as that presented in 
Fig. 2(a). In the FC case, the sample was first cooled 
from normal state (T > Tc(0)), in the presence of a cer¬ 
tain applied field {H < Hc 2 ) down to T = 2K and a 
given x'{H) value was recorded. Thereafter, the sam¬ 
ple was warmed up again to a higher T (> T'c(O)) and 
the same procedure was followed for recording another 
x'fc value in a different cooling field. A collation of 
such x'pci^) at 2K is illustrated (closed circles) in 
Fig. 5. Three different field intervals can be identified; (I) 
Across the range, 4kOe< H < the ZFC and FC 

curves can be seen to be well separated with the latter 
possessing more diamagnetic values than the former. Fol- 



Figure 5: (Color online) A comparison of x'(H) plots 
recorded in ZFC and FC modes at T = 2 K. Across the field 
range 4kOe< H < the FC state is more diamagnetic 

than the ZFC, implying the inequality, . A re¬ 
versal of this situation, i.e., is observed below 

4kOe (~ H*) in an expanded portion of x'[H) as shown in 
the inset panel. 


lowing equations (Ref. M) (i) X' ~ -1 + ahacljc and 
(ii) x' ^ —Pjc/hac, a more negative x' implies larger 
jc. Therefore, in the field interval 4kOe< H < 
the FC state exhibits a higher jc than that in the ZFC 
state (i.e., jc'" > jc^^)- As per a des crip tion of Larkin- 
Ovchinnikov collective pinning theory [30Ll3l| for weakly- 
pinned superconductors, jc relates inversely to the vol¬ 
ume (Vc) of a domain within which the vortices are col¬ 
lectively pinned and remain well correlated. Therefore, 
a higher jc amounts to a smaller 14, and hence signi¬ 
fies a more (strongly - pinned) disordered vortex matter. 
(II) The history-dependence in jc{H) tends to cease 
above as we observe the two curves to overlap 

there. This is in line with the understanding that the 
vortex matter above the peak field of the PE is gener¬ 
ally believed to be ’disordered in equilibrium’ [13 ■ (HI) 
The ZFC and FC curves appear to overlap at low fields 
(Ff < 4kOe) as well, as apparent from the two sets of 
data points in the encircled portion of x'{H) (cf. main 
panel of Fig. 5). However, a closer examination of the 
two x'{H) curves in this region on an expanded scale 
(see inset panel of Fig. 5) reveals slightly less diamag¬ 
netic values (at a given H) for the FC state than that for 
the ZFC mode, which implies the reversal of the above 
inequality; i.e., jff^ < jc^^ for H < 4kOe. We have 
marked this crossover field value (^ 4kOe) as H* in the 
inset panel of Fig. 5. Such an inequality {jff^ < jc^^) 
below H* is found to be very robust as it can be observed 
at other temperatures as well (all data not shown here 
but the H* values at different temperatures have been 
displayed in the H-T phase diagram as shown ahead in 
Fig. 8). 
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Figure 6: (Color online) Temperature - dependent x'{T) 
curves obtained in ZFC, FCC and FCW modes at Hdc = 
(a) 8kOe, (b) 7kOe, (c) 6kOe and (d) 3.5 kOe. At higher 
fields (8kOe, panel (a)), the FC curves (FCC and FCW) are 
more diamagnetic than the ZFC and hence, they imply the 
inequality . On reducing the field to 6 kOe (panel 

(c)) and then to 3.5kOe (inset panel of (d)), the ZFC curve 
becomes the most diamagnetic below a characteristic temper¬ 
ature (i.e., < jc^^)- At H — 3.5 kOe (see main panel 

of (d)) there is no signature of PE in the ZFC mode while a 
residual fingerprint of it can be seen in the FC curves. 


The H~T phase space region below H* was further 
explored via isofield x'{T) scans recorded in different 
modes, viz., zero field-cooled (ZFC), field-cooled cool¬ 
down (FCC) and field-cooled warm-up (FCW) runs as 
depicted in Fig. 6. The sample was cooled in (near) 
zero field from the normal state down to 1.8 K, and then 
a desired magnetic field was applied. Thereafter, the 
x'(T) data were obtained (marked as ZFC in Fig. 6), 
while warming the sample to higher T. Following this, 
X '(T) data were again recorded while cooling the sample 
down to 1.8 K (marked as FCC) in the presence of the 
same applied held, and thereafter while warming it to¬ 
wards the normal state (marked as FCW). At H = 8 kOe 


(Fig. 6(a)), the shielding response is found to be less dia¬ 
magnetic for the ZFC case than that observed for the 
FCC and FCW modes (i.e., j^^{8kOe) > j^^^{8kOe) 
at all temperatures). This implies that the vortex mat¬ 
ter created in the ZFC mode is better spatially ordered 
among the three modes which is consistent with the ob¬ 
servations made at this held value in Fig. 5. At a lower 
held {H = 7kOe), the ZFC curve nearly overlaps with 
the FCW; these two curves stay less diamagnetic than the 
FCC for T > 4K (cf. Fig. 6(b)). On reducing the held 
further to 6kOe (Fig. 6(c)), the FCW curve can now be 
seen to be the least diamagnetic amongst the three curves 
at all temperatures. Further, there is an unusual inter¬ 
section of ZFC and FCC curves at a certain temperature 
marked as T +, below which the ZFC is now more dia¬ 
magnetic than the FCC and FCW curves. The ZFC curve 
at a lower held of 3.5 kOe (see inset panel of Fig. 6(d)) 
can be seen to be more diamagnetic and well separated 
from the other two (FCC and FCW) curves. This would 
imply the inequality which, in turn, sug¬ 

gests that the vortex matter created in the ZFC mode 
is more disordered than that created in the FC mode. 
This situation is the reversal of that depicted in Fig. 6(a) 
and is in agreement with the conclusion drawn from the 
history - dependent x'{H) responses below H* k, 4kOe 
(cf. Fig. 5). An expanded portion of x'(r) plots at 
H = 3.5 kOe in the main panel of Fig. 6(d) shows that 
the hngerprint of PE ceases to exists in the ZFC mode, 
whereas a tiny modulation can still be noticed across 
T+ < T < Tc{H) in the case of FCC and FCW runs. 
This suggests that the vortex state created in ZFC man¬ 
ner at H = 3.5 kOe is disordered to such an extent that 
any order - disorder vortex phase transition may not be 
identihable in the temperature-dependent warm-up mea¬ 
surements. On the other hand, at the same held value 
(3.5kOe), a nascent signature a la PE observed in FCC 
and FCW modes suggests that vortex matter is compar¬ 
atively better ordered when created in the FC mode at 
3.5 kOe. 

Another demonstration of a more ordered FC state 
(than the ZFC state) at lower held values is apparent 
from the results of thermomagnetic history - dependent 
dc M-H measurements, as illustrated in Fig. 7. The 
M-H curve shown by open circles pertains to the usual 
hysteresis loop obtained in the ZFC mode at T = 2K. 
Magnetization data recorded in the FC mode {Mpc), 
i.e., after cooling the sample from a higher T (> r,(o)) 
down to 2K in the presence of different chosen held 
values have also been displayed by open triangle data 
points in Fig. 7(a). If we associate the Mpc values 
to the equilibrium (reversible) magnetization {M^q) and 
ignore the contribution due to bulk currents that can 
get set up due to a gradient in macroscopic held (H) 
(as in the case of ZFC), then, as per equation Meg = 
+ M{H~))/2 the Mpc data are supposed 
to fall in the middle of the ZFC envelope loop (as, for 
example, evident in Ref. M)- Figure 7(a) depicts an 
unexpected scenario, wherein the Mpc data points stay 















E Thermomagnetic history dependence in jc{H, T): FC state is more ordered than the ZFC at low fields 


7 


a> 

O 




a> 

O 




60 

40 

20 

0 

-20 

-40 

-60 

20 

10 

0 

-10 

-20 

-30 


1 ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 ' 

:(a) 

o ° ZFC 

A Mpc before impulse 

A Mpc after impulse 

' H* C . 

-I—I—I—|—I—I—I—I—|—I— 

0 ^ 
o o 

o o 

0 o 

o° PE °o 

O o 

f 4 * * %OJnm> 

a X — 0 ° 

▲ 1 — Z IV o 

/ Yb3Rh4Sni3 

o 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 ■ 1 1 

“y ■ 

1 1 1 1 1 1 1 1 1 1 

0 5 10 15 

20 25 

j 1 1 1 1 1 1 1 1 1 1 1 

o (b) 

I I I 1 I I I I _ 

o 

o 

T = 2K 

PE 



5 10 15 

H (kOe) 


Figure 7: (Color online) (a) M-H loop at T = 2 K in the ZFC 
mode (open circles). The field-cooled magnetization values 
{Mfc, open triangles) for 5kOe< H < 15kOe unexpectedly 
lay outside the ZFC envelope loop. An impulse of an ac field 
(amplitude = 10 Oe, / = 211 Hz) imposed on the FC states 
created at different H (open triangles) leads the Mfc values 
to fall well within the loop (closed triangles), (b) Minor hys¬ 
teresis loops at T = 2 K traced while reducing the field to 
zero after creating the FC states at JT = 5, 10, 15 kOe and 
imposing an ac field impulse on them. 


outside the envelope hysteresis loop for field (H) values 
in the interval, 5kOe< H < 15kOe. We can try to 
rationalize this using a model due to Clem and Hao [s^ 
which accounts for the FC magnetization of a type - II su¬ 
perconductor when a gradient in macroscopic field gets 
established as a consequence of flux expulsion during the 
FC mode. In such a situation, the Mpc values deviate 
from Meq (ideal case of no flux - pinning). This deviation 
of FC magnetization from Meq values is indeed antici¬ 
pated to be governed by the strength of bulk pinning 
at a given field; the extent of deviation is more (less) 
for stronger (weaker) pinning strength (jc) (see Fig. 5 in 
Ref. [1^). The Mpc values (see open triangles, Fig. 7(a)) 
falling outside the envelope loop in the present case, for 
H > 5 kOe are due to a larger jc value associated with 
the vortex matter created in the FC mode and the flux 
density gradient [s^ that gets set-up while field - cooling. 
Note that below H = 5kOe (close to H* ~ 4.5kOe), the 


FC state seems to have a lower jc value (better ordered) 
as the Mpc value here remains well within the envelope 
(i.e., nearer to Meq). 

As a further experimentation, on each occasion, af¬ 
ter a FC state was created at a chosen field value, we 
momentarily perturbed it with an ac field impulse of 
amplitude 10 Oe at a frequency of 211 Hz (applied for 
about 6seconds). It is curious to note that the mag¬ 
netization recorded after the impulse treatment (closed 
triangles data points in Fig. 7(a)) fall well within the en¬ 
velope loop. It appears as if the magnetization response 
of a FC state after perturbation by an impulse conforms 
to the anticipated respective equilibrium value (Meq) 
(closed triangles data points) located in the middle of 
the ZFC envelope loop. The imposition of an ac field im¬ 
pulse results in the reconfiguration of an unperturbed FC 
state into state, comparison of whose jc value with that 
of the corresponding ZFC state can be very instructive. 
We thus recorded the ’FC minor hysteresis loops’. Fig¬ 
ure 7(b) shows the magnetization curve obtained while 
reducing the dc field to zero after the (perturbed) FC 
states were created (encircled) at H = 5kOe, lOkOe 
and I5kOe. The initial magnetization value recorded 
while decreasing the field from H = 5kOe undershoots 
the ZFC envelope, and the magnetization curve there¬ 
after traverses a path that remains within the envelope 
loop as the field is ramped down to the zero value. Tak¬ 
ing cue from the linear relation between hysteresis width 
and jc [13 , the (perturbed) FC state is reckoned to have 
a lower jc value than that in the ZFC for fields below 
H*. This observation fortifies the inferences drawn from 
data in Figs. 5 and 6. At H = lOkOe, the magnetiza¬ 
tion data, while reducing the field, first overshoots the 
ZFC envelope though marginally, and thereafter it re¬ 
traces the reverse leg (M (H~)) of the ZFC magnetization 
curve. It can be argued that the FC state here exhibits jc 
value slightly higher than that in the corresponding ZFC 
states. A significant overshooting of the ZFC envelope 
loop is witnessed at a higher field of 15 kOe, which indi¬ 
cates a larger jc value for the FC state than that created 
in the ZFC mode. Note that the field value of 15 kOe lies 
in the anomalous region, iJ°"'’“^(r) < H < 
of the H-T space at 2K (Fig. 4) which, as discussed in 
section-IV ahead, is ’partially disordered in equilibrium’ 
as the multi-domain VG state. 


IV. DISCUSSION 

The magnetization measurement techniques (both ac 
and dc) employed in the present work have lead to new 
revelations in Yb 3 Rh 4 Sni 3 which have been summarized 
in the form of a modified H-T phase diagram in Fig. 8. 
These include, (i) the occurrence of a broad anomaly lo¬ 
cated deeper in the mixed state triggering at i7°”’“‘’(r) 
line as seen in the x'{H) curves (Fig. 2), (ii) identifi¬ 
cation of onset position of the PE at 77°"’'*'”, line 

as obtained from both the dc M-H (Fig. 1) and the ac 
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x'{T) (Fig. 3) plots, and (iii) a characteristic line H*{T) 
extracted from Fig. 5, below which a vortex matter cre¬ 
ated in the FC mode is somewhat more ordered than 
that obtained in the ZFC manner. Following a descrip¬ 
tion of L - O theory HO,mi, we argue that the anomaly in 
x'{H) (Fig. 2) reflects a shrinkage in 14 as jc oc IjyjVc 
and hence, can be termed as an order-disorder transi¬ 
tion in the vortex matter. Note that the L-0 theory 
is well applicable for weakly-pinned superconductors as 
in the present case of Yb 3 Rh 4 Sni 3 (ratio of depinning 
and depairing current densities is of the order of 10“^). 
The location of i?°"’“°(T) line deeper in the mixed state 
and its weak temperature - dependence outside the boxed 
region in Fig. 4 can be argued to relate to the disorder- 
induced transition [l|, 0,0,0] of a quasi-ordered (elastic) 
vortex lattice a la BG phase to the dislocation-mediated 
(multi - domain) VG phase. Note that the peak field of 
the anomalous variation in jdH) as reflected in dc M-H 
and ac x'{H) plots, viz., ■'^'^(T) (closed squares) and 
}jpk,ac^rp^ (open squares), is consistent, as the two sets 
of data points extracted from two different measurements 
can be seen to be overlapping in Fig. 8. This would im¬ 
ply that the complete amorphization of the vortex mat¬ 
ter (which is generally believed to occur at peak field of 
the PE [U) occurs at the same field/temperature value 
while performing measurements via ac and dc magneti¬ 
zation techniques. 

It has been shown earlier that the annealing effects 
on the vortex matter, produced either by an ac driv¬ 
ing force llTl I 33 I or by repeated cycling of the dc mag¬ 
netic field [ 34 L l35l| , eliminate either partly or completely 
a (metastable) disordered vortex phase yielding an or¬ 
dered vortex configuration. Also, the residual presence 
of a disordered (metastable) vortex phase usually governs 
the onset position of the order-disorder (PE) transition(s) 
[^[^,1^ . As a usual behavior, the onset position of the 
PE shifts to higher field values after the annealing of (dis¬ 
ordered) vortex phase, whereas the same moves to lower 
fields when the amount of quenched disorder is larger 
[U [ 2 ^ . Surprisingly, the shaking effect of an ac driving 
force appears to be counterintuitive in the present study, 
as it promotes the spatial disordering of the vortex matter 
rather than its usual role of improving the state of spatial 
order. This is well apparent by the (early) occurrence of 
an order - disorder transition at FJ°"'’“'^(T) line (cf. Fig. 
8) in the x'(iF) runs which involve the ac drive. In the 
absence of an ac driving force as in dc M-H scans, the 
onset of order-disorder transition (PE) is seen at higher 
fields; iF°"’‘^'^(T) values in Fig. 8 are located significantly 
above the BG to VG transition (iF°"’“'^(T)) line. If one 
assumes that the shaking effect on the vortex array by 
an ac drive results in the lowest equilibrium-like state of 
the system (as in Ref. [13) under given circumstances, 
then the region bounded by the phase lines, iJ°"’’“'^(T) 
and in Fig. 8 can be accepted as “disordered 

in equilibrium ” in the form of multi-domain vortex glass 
phase. The present hndings echo similar assertions made 
recently by us [l^ in another study in a single crystal of 



Figure 8: (Color online) Complete vortex phase diagram in 
our crystal of Yb 3 Rh 4 Sni 3 . Various regions have been identi¬ 
fied (refer text for details). A delineation between the BG to 
VG transition (occurring at //“"’““(T)) and the PE anomaly 
(at Hp"’‘^“(T)), is well apparent. The peak field values 
qjp^^d.c^T'j and obtained respectively from the M- 

H and x'{H) data almost fall on the same phase boundary. 
The region between /Hf'“dT)) and the HirdT) 

line presumably comprises pinned amorphous matter while 
the narrow space between Hirr{T) and Hc 2 {T) may involve 
an unpinned amorphous vortex matter. Below H*(T) (ob¬ 
tained from Fig. 5), the H-T phase space region comprises 
(reentrant) disordered vortex matter, a vortex configuration 
created here in the FG manner is found to be slightly more 
ordered than that in the ZFG mode. 


a low Tc superconductor, Ca 3 lr 4 Sni 3 . However, in that 
study in Ca 3 lr 4 Sni 3 , a distinct demarcation between the 
BG to VG transition line and the locus of the onset of the 
PE anomaly in its vortex phase diagram was not appar¬ 
ent, which has now been clearly sorted out in the phase 
diagram of Yb 3 Rh 4 Sni 3 . 

The nature of BG to VG transition had been argued 
to be of first - order in a high Tc superconductor. To 
fortify this proposition, we may emphasize that the re¬ 
gion bounded between and in Fig. 8 

which is shown from the y'(iF) data to be disordered 
in equilibrium is indeed an ordered vortex phase when 
viewed from the outcomes of the dc M-H measurements. 
Therefore, vortex phase in this region can be treated as 
a superheated ordered BG phase. This attests the first- 
order nature of the BG to VG transition in a low Tc 
superconductor. 

The change in magnetic field tunes the inter-vortex 
spacing which, in turn, can influence the balance between 
the strength of vortex pinning and the (elastic) interac- 
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Figure 9: (Color online) \x'lx'v\ data (normalized with re¬ 
spect to X' value at peak position of the PE) plotted against 
the reduced temperature, T/Tp at various constant dc fields. 
The shape of the PE anomaly changes with field. At // = 
4.5 kOe, there is a tiny peak feature while a much prominent 
peak feature can be noticed at a higher field, H = 7kOe. 
The PE feature gets broadened at a further elevated field of 
15 kOe. 


tions between the vortices. As a result, the extent of spa¬ 
tial ordering/disordering in the vortex matter may vary 
in different regions of the H-T phase space. Further, an 
enhancement in effective pinning at a given field value 
can give rise to a qualitative chan g^e i n the evolution of 
the size of the Larkin domain , as had been ex¬ 

perimentally demonstrated via the evolution of PE 
feature in the isofield x'{T) measurements in crystal(s) 
of 2 H-NbSe 2 . Motivated by this, we shall now exam¬ 
ine how the PE feature evolves with magnetic field in the 
case of our crystal of Yb 3 Rh 4 Sni 3 . We show in Fig. 9 the 
normalized |x VXpl plots (x'p corresponds to x' value at 
peak position of PE where vortex matter is most disor¬ 
dered for a given inter-vortex spacing and the underlying 
quenched random disorder) against the reduced temper¬ 
ature = T/Tp, [Tp is peak temperature of PE obtained 
in the ZFC mode for various fixed fields). The PE fea¬ 
ture remains absent at fields below H = 4.5 kOe (data 
not shown here), which may be because of a disordered 
vortex matter prevailing there. Note that the larger inter¬ 
vortex spacings at low flux density results in weaker in¬ 
teraction between the vortices to counter disordering in¬ 
fluence of the pinning centres. The vortex configuration 
created at low fields ought to be highly disordered akin 
to the reentrant a morp hous vortex matter anticipated in 
earlier studies [39l - [4^ . With the increase in field, the 
interaction effects strengthen and consequently, a sparse 
disordered vortex configuration can progressively change 
into an ordered one, as has been demonstrated earlier in 
Refs. [il,|U. This is apparent from Fig. 9 by the appear¬ 
ance of a tiny PE feature firstly at H = 4.5 kOe. A more 
pronounced signature of PE occurring at H = 5.5 kOe 


suggests further improvement in the spatial ordering of 
vortex matter with the increase in field prior to the on¬ 
set of PE. At a higher field of 7kOe, there emerges a 
very well developed PE feature, which is indicative of a 
well-ordered vortex (BG) phase prevailing prior to the 
onset of the PE transition at this field value. A dif¬ 
ferent scenario is, however, depicted at a higher field, 
H = 15kOe, wherein the order-disorder transition can 
be seen to be much broader. We surmise that the vor¬ 
tex phase prior to the order-disorder transition at this 
field value is partially disordered (a la multi-domain VG 
phase) so that any order-disorder transition would reflect 
a lesser shrinkage in 14 (as compared to that at lower 
fields) or equivalently, lesser increment in jc dsim as¬ 
sociated with this state. In all, the \x'/x'p\ plots of Fig. 
9 have lead us to infer that a disordered vortex phase 
at low field end {H < 4.5 kOe) transforms into an or¬ 
dered one (BG phase) with increase in field. A further 
enhancement in field may transform an ordered vortex 
phase (BG) into a partially disordered (VG) state prior 
to the onset of PE. 

It is useful to recall now a recent study [i^ employ¬ 
ing the small angle neutron scattering (SANS) technique 
in another crystal of Yb 3 Rh 4 Sni 3 (which would typi¬ 
cally have different amount of quenched random disor¬ 
der) providing evidences at microscopic level in support 
of our experimental results. In the SANS study, the 
evolution of spatial order in vortex matter with vortex 
lattice spacing (oq) has been well-depicted by plots of 
field-dependences of the normalized correlation lengths 
(both longitudinal, ^l/oq and transversal, ^ r/a n) as well 
as their ratio (Cl/Ct) (see Fig. 3 in Ref. [4^). These 
lengths relate to the longitudinal and radial dimensions 
of the volume within which the flux lines remain corre¬ 
lated. The full width at half maximum of the rocking 
curves in longitudinal and transverse directions reflect¬ 
ing the apparent quality of the spatial order and that of 
the Bragg peaks remains nearly flat in the field interval 
hkOe < H < 16 kOe (cf. Fig. 2 in Ref. [1^). The 
normalized transversal correlation length (Ct/oq) starts 
decreasing more rapidly as the field reduces below 4 kOe 
[ 4 ^. Near 4.5 kOe, where we just start to observe the PE 
feature in Ix'iT)\ scans in our sample, the said ratio has a 
value ^ 2 . This could imply that when the quality of 

spatial order in the vortex matter gets more compromised 
(i.e., Ct/oq < 2), the PE feature would not surface up. 
The longitudinal correlation length was seen [4^ to be 
nearly fifty times longer than It has been argued 
that the ratio /^t reveals field-induced changes in the 
parameter a/ cm/cqs, where C44 and cee are respectively, 
the tilt and shear modulii of hexagonal vortex lattice. 
This ratio in their sample remains nearly constant be¬ 
tween 5kOe and 12.5 kOe and starts to fall above it and 
the rate of decrease becomes very steep above 17kOe (see 
inset of Fig. 3(a) of Ref. M)- Such a trend could be 
interpreted that the quality of spatial order undergoes 
a change above about 12.5 kOe and the elasticity of the 
vortex lattice plummets above 17.5 kOe. In the context of 
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our sample, this could rationalize that BG to VG transi¬ 
tion happens above 12.5 kOe (see Fig. 8). The inset panel 
in Fig. 3(b) of Ref. shows that the simulated (nor¬ 
malized) vortex-vortex interaction in their sample does 
not continue to monotonically rise with magnetic field 
as the inter-vortex spacing progressively decreases with 
increase in field. The vortex-vortex interaction (for the 
Yb 3 Rh 4 Sni 3 crystal investigated in Ref. El) reaches a 
maximum value near 12kOe and starts to decrease above 
about 14 kOe thereby implying that the onset of this de¬ 
crease is a consequence of notion of enhancement in ef¬ 
fective pinning at large field values. Near 20kOe, the 
value of the normalized vortex-vortex interaction force 
for this compound is same as at about 6kOe [i^. The 
rapid decrease in this parameter below 5 kOe reflects the 
dominance of pinning in the field domain of reentrant 
amorphous state (cf. Fig. 8). 

In our present study, the \x'lx'p\ plots of Fig. 9 sum¬ 
marize the evolution in the spatial order of vortex mat¬ 
ter with magnetic field, which ties up well with the re¬ 
sults of SANS study [4^. The extent of spatial order¬ 
ing/disordering in the vortex matter at a given field can 
be quantified by the ratio of |x VXpl value just prior to 
the PE and that at the peak position (see Tp in Fig. 9) 
of the PE. Note that the correlation length ought to be 
minimum at the peak (Tp) of the PE and maximum (for a 
well ordered vortex lattice) just before the PE. In Fig. 9, 
the said ratio of \x'lx'p\ at onset and peak positions of 
the PE can be seen to be increasing smoothly with in¬ 
creasing magnetic field for 4.5 kOe < H < 7 kOe, which 
is consistent with the variation of correlation lengths with 
field for H < lOkOe, as anticipated in Ref. [^. This 
would imply that a very well-ordered BG phase emerges 
nearly at H = 7kOe from a disordered state prevailing 
below 4.5 kOe. Further, at an elevated field (15 kOe), this 
ratio decreases substantially as if a transition from the 
BG phase into a VG like phase has occurred there. 

In the end, we draw attention towards the thermomag- 
netic history effects investigated at lower fields (Figs. 5- 
7). It is customary to witness a more disordered FC 
state than the ZFG in the vicinity of the PE phenomenon 
prl - l^ . However, the history-dependent magnetization 
behavior at lower fields {H < H*) studied here in 
Yb 3 Rh 4 Sni 3 have revealed that the vortex matter in the 
ZFG mode is more disordered than that in the FG case. 
The following scenario seems plausible to explain this fea¬ 
ture. In the ZFG mode, the injection of vortices at high 
velocities into a superconducting specimen is influenced 
by surface barriers and edge effects [H, such that 
the vortices eventually penetrate the sample through the 
weakest point of the barrier. Inside the superconduct¬ 
ing specimen, the vortices moving in bundles get ran¬ 
domly pinned at respective pinning sites such that the 
inter-vortex spacing is non-unique which leads to a non- 
uniform distribution of the flux-density. As mentioned 
earlier, the vortices are well separated at lower fields and 
therefore, the (elastic) interactions between them remain 
weak resulting in a stronger pinned (disordered) vortex 
configuration in the ZFG mode. Although the vortex 


matter created in the FG mode also remains disordered 
at lower fields, however, the vortices in this mode nucle¬ 
ate more uniformly as they remain oblivious to disorder¬ 
ing effects of the surface barriers and the non-uniform 
injection through edges, etc. Therefore, at lower fields, 
the vortex configuration in FG mode is less disordered 
than in the case of ZFG. At higher fields, as the vortices 
get closer to each other and hence, the interaction effects 
become dominant yielding a better spatial ordering even 
during the moving state of creation of vortex matter in 
the ZFG mode. On the other hand, during the FC mode, 
the vortex matter has to cross the PE boundary in the 
H-T phase space which results in supercooling the dis¬ 
ordered vortex matter prevailing at the peak position of 
the PE anomaly. Thus, at higher fields, the vortex mat¬ 
ter created in the ZFG mode is more spatially ordered 
than that in the FC mode. 


V. CONCLUSION 

We have investigated via magnetization measurements, 
a weakly - pinned single crystal of a low superconduc¬ 
tor, Yb 3 Rh 4 Sni 3 . The present results have lead to the 
identification of BG to VG transition line and sketch of a 
characteristic field {H*{T)) line in the H-T phase space 
of this compound. Surprisingly, the SMP transition like 
phase boundary has been unearthed under the combined 
influence of an ac driving force and a continuous magnetic 
field sweeping involved in the (isothermal) ac susceptibil- 
ity (x '(77)) measurements. This transition was, however, 
not observed in the (isothermal) dc M-H loops and the 
temperature - dependent ac susceptibility scans (x'(T)). 
The latter two modes of measurements yield signature 
of only the quintessential PE anomaly signaling the col¬ 
lapse of elasticity of vortex solid at higher fields. An 
apparent demarcation between the domain of SMP like 
transition and the onset of PE anomaly has been made 
in the vortex phase diagram of Yb 3 Rh 4 Sni 3 . The results 
presented in our specimen answer in affirmative the ques¬ 
tion of generic nature of BG to VG transition in a pinned 
superconductor. These also find supports from another 
experimental study [45l| at microscopic level in a sample 
of the same compound. In the low field region {H < 77*), 
the vortex matter is construed to be highly disordered. 
Here, a vortex state created in the FC mode is found to 
be more ordered than that obtained in the ZFG mode. 
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